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Reliability  in  a  Random  Environment 

Edward  B.  Rockower 
Naval  Postgraduate  School 
Dept,  of  Operations  Research 
Monterey,  CA.  93943 

ABSTRACT 

Correlations  in  the  failures  of  subsystems  or  components  may  arise 
when  they  share  a  common,  random  environment.  We  show  that  a  natural 
tool  for  calculating  the  reliability  of  such  systems  is  the  characteristic 
functional  of  the  random  hazard  rate.  h(t).  Some  general  results  for  the 
reliability  of  series  and  parallel  systems  in  terms  of  the  characteristic 
functional  of  the  hazard  rate  are  derived  and  applied  to  a  number  of 
models  of  random  environments.  The  applications  include  random  hazard 
functions  arising  from  I)  non-fatal  shocks  of  random  amplitudes,  2)  a 
Markovian.  Gamma-marginal  stochastic  process,  3)  system  stress  related 
to  un-repaired  damage  from  incoming  rounds,  and  4)  impulsive  and 
accumulated  heat  stresses  from  a  laser  battle. 

I  INTRODUCTION 

Consider  a  system  whose  components  are  subject  to  a  random 
environment1,  possibly  including  random  shocks.  In  general,  the  ith 

component  hazard  r3te,  h,(t),  may  be  written  as  the  sum  of  tw:  parts. 

The  rirst  (deterministic)  part.  X,(t).  accounts  for  wearout  and  those 
random  fa1  lures  which  occur  independently  for  each  component.  The 
second  (stochastic)  part,  g,(tj,  arises  from  those  environmental 
conditions  which  are  shared  by  two  or  more  components  in  the  system. 


Aocessloo  For 

mis  ORAM 

DTIC  TAB 

Unannounced 

% 

□ 

Justification _ 

By 

Distribution/ 


Availability  Codes 
jAvail  and/or 


Dist 


M 


Special 


We  develop  a  general  formalism  for  calculating  the  reliability  of  a 
system  in  a  random  environment  in  terms  of  the  characteristic 
functional  of  the  hazard  rate,  various  models  of  random  environments 
are  then  proposed  and  the  corresponding  expressions  for  the 
characteristic  functionals  are  calculated.  The  reliability  of  systems  in 
such  environments  are  then  found  using  these  characteristic  functionals. 
The  equivalent,  deterministic  hazard  functions  for  single  components  in 
some  of  these  models  are  found  not  to  be  constant  in  time2. 

In  Section  II  the  general  formalism  is  introduced  and  illustrated  for  a 
single  component  subject  to  an  environment  which  includes  non-fatal 
shocks  with  random  amplitudes,  in  Sections  III  and  IV  the  reliability  of 
series  and  parallel  systems  is  presented  in  general  terms  and  applied  to 
the  non-fatal  shock  model.  Section  V  presents  a  model  for  the  hazard 
rate  as  a  Markovian.  Gamma-marginal  stochastic  process.  Finally,  in 
Section  VI  we  present  three  models  in  which  the  failure  rate  is 
proportional  to  system  stress.  Each  of  these  models  can  be  applied  to 
multi-component  systems  using  the  methods  developed  in  the  earlier 
Sections. 


One  Component 

The  simplest  case  is  that  of  one  component  with  hazard  rate  h(t),  a 
non -negative  stochastic  process.  The  reliability  of  the  component, 
subject  to  a  particular  realization  of  the  environment  (hence,  h(t)  )  is 
def  med  as. 

anit)  =  Pr(T>t]n.  (1) 
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where  T  is  the  random  variable  (r.v.)  equal  to  the  time  the  component 
fails  and  the  subscript  denotes  a  particular  realization  of  the  stochastic 
process,  h.  Using  the  usual  definition  of  the  hazard  rate3,  this  is. 

t 

Rh  =  exp[  -Jh(t‘)dt‘  ].  (2) 

o 

Hence,  averaging  over  the  random  environments  we  have, 

t 

R(t)  =  Eh[  Rh  1  =  E{  exp[  -Jh(t‘)dt*  ] }  (3)' 

o 

or. 

R(t)  =  Ch[Tl(-)]  l^j  ,  (A) 

where, 

t 

Cn  t[T|(-)]  =  E{  exp[  i jTj(t’)  h(t')dt'  ] },  (5) 

o 

defines  the  characteristic  functional4,5  of  the  hazard  rate.  In  other 
words,  this  characteristic  functional,  evaluated  for  a  particular  value  of 
the  test  function,  T}(t),  yields  directly  the  reliability  of  a  component. 

We  will  see  below  how  this  generalizes  to  multi-component  systems  in  a 
random  environment. 


We  can  write  the  reliability  of  a  component  or  system  in  terms  of  an 
effective,  deterministic  hazard  rate,  h(t),  using  the  usual  definition  of 


hazard  rate. 


h(t)  =  -  d/dt  in  {  R(t) }, 

but  we  must  remember  not  to  use  this  “hazard  rate”  to  calculate 
reliability  of  larger  aggregations,  unless  such  aggregations  are  composed 
of  components  or  subsystems  not  sharing  a  common,  random  environment. 

Deterministic  hazard  rate  plus  random  shocks 

If  the  environment  gives  rise  to  non-fatai  random  shocks  with  rate  tf(t), 
then  the  hazard  rate  may  be  represented  by  a  deterministic  function, 

X(t).  plus  delta  function  spikes,  possibly  with  random  amplitudes.  For 
convenience  we  will  sometimes  replace  X(t)  or  £(t)  with  constants. 
However,  in  all  the  reliability  formulas  derived  here,  Xt  can  always  be 
replaced  with  JX(t’)dt’  for  non-constant  background  hazard  rate,  and  3t 
can  be  replaced  with  Jff(t')dt'  for  a  non-homogeneous  shock  process. 

Note  that  X(t)  may  be  really  an  effective  hazard  rate,  obtained  from  the 
preceding  equation,  when  one  portion  of  the  random  hazard  rate  is  unique 
to  each  individual  component,  i.e.  not  shared  among  separate  components. 

The  hazard  rate  for  non-fatal  shocks  occurring  at  times  tj  is, 
n 

h(t')  =  X  ♦  2  o<j  S(t*  -  tj).  (6) 

]  =  ’ 


A 
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where  the  pdf  of  the  tj  is  2f(t)/n  ,  n  =  |2f(t’)dt\  and  n  is  a 

o 


Poisson  r.v.  with  mean  n  . 

The  characteristic  functional  of  this  process,  derived  in  Appendix  A,  is, 

t  t 

Chth(-)1  5  exp{  »*•  J*n(t*)dt*  ♦  J 2f(t’)  [Cj  Tj(t’) )  -  1  ]dt  ).  (7) 

o  o 

where  C(>  is  the  ordinary  characteristic  function  of  the  shock 

amplitudes,  <*j.  In  particular,  if  the  shock  amplitudes  have  an 
exponential  distribution  (parameter  S)  then, 


C^Coj)  -  S/[S  -  i  ujI. 


(8) 


and. 


t  t 

Ch  t[Tj(  ■)]  =  exp{  iX  /T](t')dt'  ♦  iftf(t')  Tj(t')/[6  -  iTj(t')]  dt'l  (9) 

o  o 

Hence,  setting  t\(1)  =  i,  the  reliability  is, 

R(t)  =  exp{  -Xt  -  n  /IS  ♦  1] },  00) 

or,  assuming  a  stationary  environment  (2f  =  const.)  and  using  c<  -  1/8, 


R(t)  -  exp(  -Xt  -  tft  cx  /[cx  *  I] }. 


(ID 


The  effective  hazard  rate,  h(t),  is  constant  for  this  model.  Note  that  as 
o<  ->  oo  the  rate  of  failure  becomes  X  plus  the  rate,  tf,  of  occurrence  of 
(fatal)  shocks  (in  fact,  X  could  be  a  background  rate  of  infinitely  high 
shocks). 

111.  Two  different  components  in  the  same  environment 
If  two  or  more  components  are  in  different,  independent  environments 
the  system  reliability  follows  from  the  usual  formulas3  expressed  in 
terms  of  the  individual  reliabilities.  We  consider  here  two  different 
components  seeing  the  same  environment.  The  reliability  of  two 
components  in  series  subject  to  a  given  realization  of  the  environmental 
conditions  is. 

Rh(t)  =  Ph[T>t)  =  Rh|(t)  •  Rh2(t).  (12) 

If  the  components  are  in  a  parallel,  redundant  system  the  reliability  is, 

Rh(t)  =  1  -[I  -Rhl(t)l-[!  -Rh2(t)l  (13) 

-  exp[  -Jh^Odt'  ]  ♦  exp(  -Jh2(t’)dt’  1  -  exp{  -Jlh^t’)  ♦  h2(t‘)l  dt’ }, 

where  all  integrals  are  0  ->  t.  Hence,  averaging  over  realizations  of  the 
environment  wa  have,  using  the  definition  of  the  characteristic 
functional. 

R(t)  =  (14) 

ch,t^1'^  U(  )='  +  ch,tW'^  I  t\(-)=i  "  clvh2,tW'^ 


Clearly,  the  last  term  also  represents  the  reliability  of  a  series  circuit 


with  the  same  two  components.  In  the  following  it  will  be  understood 
that  all  characteristic  functionals  are  for  the  processes  over  the 
interval  [0,  t]. 

As  an  example,  consider  two  components  in  series,  seeing  the  same 
shocks,  but  experiencing  different  amplitudes,  c*,  k£;  where  <x  and  £  are 
i.i. d.  exponential  r.v’s.  Then  h,  ♦  h2  is  given  by 

n 

h(t')  =  X,  *  X2  ♦  I(cxj  ♦  k$j)8(t' -  tj).  (15) 

j=1 

In  this  case  in  Eq.  (7)  is  replaced  by  the  characteristic  function  of 
the  sum  of  two  exponentials, 

Ca(a>)  -  8/lS  -  i(i)l  •  S/[S  -  iu)k),  (16) 

and  X  is  replaced  by  X,  +  X2,  yielding, 

t 

Ch,*h2[Tl(')l  =  exp{  i(X1+X2)J,n(t‘)dt‘ }  •  (17) 

o 

t 

exp{  J2f(t')dt'[  8/{S  -  iTj(t')]  •  8/[8  -  iT](t  )ki  -  i)  }. 
o 

The  reliability  of  the  series  system  is  then  found  by  substituting 
Tj(t')  =  i ,  yielding. 

Rseries(t)  =  exP{  '(VM  -  *t[  l  -  82/[(S  ♦  0(8  ♦  k)1 1 ).  08) 

The  parallel  system  reliability  is  given  by, 


R||(t)  -  R|(t)  +  R2(t)  -  ^series' 


(19) 
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(t). 

where  R ]  2  are  the  same  as  Eq.  (10)  with  8  replaced  with  S  and  8/k, 
respectively. 

if  cx  =  3  in  the  above,  i.e.  the  amplitude  of  shocks  seen  by  the  two 
components  are  proportional,  then  we  have, 

n 

h,(t‘)  +  h2(t‘)  =  X,  ♦  X2  ♦  (l+k)  X  <xj  5(t*  -  tj).  (20) 

j=> 

This  is  the  same  as  for  a  single  component,  except  X0  ->  X}  +  X2  and 
8  ->  8/(l*k).  hence. 

t 

chl+h,h(-)]  =  exp{  i(X;*X2)/r|(t')dt' }  •  (21) 

t  0 

exp{  JWJdt’I  S/[S  -  iT|(f)  (k*l)]  -111. 

0 

and  the  reliability  of  a  series  circuit  is, 

Rsenes(t)  =  exp{  -  (X,  +  X2)t  -  2ft  [(k+ 1  )/(8+k+ 1 )]  }.  (22) 


The  reliability  of  a  parallel  circuit  is  again  given  by  Eq.  (19). 


Rh(t)  =  Rh,^  V°  V°‘  " 

t 

=  exp{  -J(h,  ♦  ■  •  •  ♦  hm)df  }.  (23) 

o 

Hence,  averaging  over  the  environment,  the  reliability  is, 

s  w -  (24) 

a)  If  all  components  respond  to  a  given  shock  with  independent 
amplitudes  we  obtain, 

Ra(t)  =  exp{  -mX0t  *  2ft  [Sm/(S+1)m  -  1]  },  (25) 


b)  if  the  hazard  rate  is  exactly  the  same  for  all  components, 

Rjj(t)  =  exp{ -m\0t  +  2ft  [6/(S+m)  -  I]}.  (26) 

The  effective  hazard  rate,  h(t),  is  again  constant.  Note  that 
(M/S)m  =  1  ♦  m/6  *  •  • 

hence, 

(M/8)m  >  1  ♦  m/S, 
or, 

£m/(S+ 1  )m  <  S/(S+m), 

which  implies  that  Ra  <  Rb.  In  other  words,  as  expected,  the  reliability 

of  a  series  system  is  higher  in  the  more  highly  correlated  environment. 
Also,  using  Eg.  (19).  it  is  clear  that  the  reliability  of  a  parallel  system 
will  oe  lower  in  a  more  highly  correlated  environment. 


The  reliability  of  m-identical  components  in  a  parallel  redundant  system 
subject  to  a  given  realization  of  the  environment  is. 

m  t 

Rh(t)  -  I  -  TT  [  1  -  exp{-Jhj(t')dt'}  ].  (27) 

j=l  o 

This  can  be  expanded  using  binomial  coefficients,  Cm^  ,  as, 
m 

Rh(t)  =  l  -!(-))  cm,  !'Rhl(t)~Rhj(t).  (28) 

1=0 

Hence,  averaging  over  the  environment,  we  have, 
m 

RCt)  =  !  -  I  (-))  e'Jxot  exp{  2ft  [si/(S*l)J  -  I] },  (29) 

j=0 

for  independent  response  to  the  shocks,  and, 
m 

R(t)  -  1  *  I(-)j  Cmj  e‘jxot  exp{  2ft  [6/(6*  i)  -  1]  },  (30) 

j=0 

if  there  is  exactly  the  same  hazard  rate  for  all  components  (each 
component  sees  exactly  the  same  amplitude  shocks). 

In  general,  for  m  identical  components  in  parallel,  all  with  the  same 
hazard  rate,  h(t)  =  X0(t)  +  g(t)  (where  g(t)  is  random), 


m 


R(t)  =  1  -  K-)i  Cmj  .  exp{ -jJX0(f)dt‘}  ■  Cg[ri(  ))| Tl(.)=j.,  (31) 

j=0 

From  the  foregoing  it  is  clear  how  to  generalize  the  formulas  for 
reliability  in  a  random  environment  to  more  general  configurations 
containing  components  both  in  series  and  in  parallel: 

1)  using  the  usual  rules  for  probabilities,  write  the  reliability  of  the 
system  in  terms  of  the  individual  component  hazard  processes, 
conditioned  on  a  given  realization  of  the  environment,  hence  of  the 

fh,(t)}  (cf.  Eq.  27). 

) 

2)  average  over  the  environment,  hence  over  the  {h j(t)}. 

3)  express  the  resuit  in  terms  of  the  appropriate  characteristic 
functionals  (cf.  Eq.  31),  and  finally, 

4)  obtain  the  effective  component  or  system  hazard  rate,  ti(t). 
v.  Exponent  i  a  I/Gamma  Hazard  rate 

Consider  an  environment  giving  rise  to  a  Markov  hazard  rate  process  with 
a  Gamma  marginal1,6,7  distribution  which  is  common  to  all  components 
.n  the  system  of  interest.  We  show  in  Appendix  B  that  the 
characteristic  functional  for  such  a  process  is, 

t  t  t" 

CghG)!  =  exp{ufdt' <f(t'i[S/lS  -  i  fdt"'?t(t")exp(-/2)'(r)dtr)]  -  I  1).  (32) 
o  t'  t’ 
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The  correlation  function  for  this  process  in  the  stationary  case 
(d  const.)  is  p  =  e"^r.  When  u=i  the  process  is  exponential. 


The  reliability  of  a  single  component  with  hazard  rate  h(t)  =  X(t)  +  g(t) 
is  then, 

P(t)  =  exp{ -JX(t*)df)  ■  Cg['n(  )||n(.)=|  .  (33) 

carrying  out  the  resulting  integrations  yields, 

R(t)  =  (34) 

exp{  -|X(t')df  -otft/(2S*l)}  • 

The  effective  hazard  rate,  h(t),  is  clearly  not  constant  for  this  model. 
This  can  be  generalized  to  multi-component  systems  with  the  methods 
from  the  previous  sections.  For  example,  the  reliability  of  m  identical 
components  in  series  is  obtained  by  replacing  X  ->  mX  and  using 
T}(-)  =  m  i  in  the  characteristic  functional.  From  Eq.  (32)  it  is  easily 
seen  that  this  leads  to  replacing  S  with  S/m  in  Eq.  (34). 


VL  Failure  Rate  Proportional  to  System  Stress 
In  this  section  we  will  consider  three  related  models  of  hazard  rate 
processes,  in  the  first  model  we  assume  the  rate  of  failure  is 
proportional  to  the  number  of  customers  using  the  system,  for  example 
the  rate  of  wear  on  a  highway  may  have  a  component  which  is 
proportional  to  the  number  of  automobiles,  N(t),  using  the  highway.  In 
the  second  mode!  the  failure  rate  of  a  major  system  (for  example  on  a 
battleship  during  combat)  is  increased  proportional  to  the  amount  of 
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un-repaired  damage  (the  number  of  hits  not  repaired  is  N(t).  the  amount 
of  damage/hit  is  some  positive  r.v.).  Both  of  these  models  are  related 
to  a  generalized  M/G/oo  queue.  Finally,  we  model  the  reliability  of 
electronic  systems  in  a  laser  battle  scenario.  Failures  may  be  caused  by 
impulsive  stress  caused  by  a  laser  hit.  or  by  accumulated  heat  from  the 
laser  hits.  It  is  found  that  none  of  the  models  in  this  section  lead  to  a 
constant  effective  hazard  rate. 

Highway  model 

Consider  a  large  multi-lane  highway  in  which  the  traffic  level  is  such 
that  cars  do  not  interfere  with  each  other.  The  number  of  cars  on  the 
highway  can  be  modeled  as  an  n/G/oo  queue.  The  rate  of  arrivals  of  cars 
to  the  highway  is  2f(t)  and  the  pdf  of  the  time  spent  on  the  highway 
section  of  interest  is  b(r)  (the  "service”  distribution).  The  hazard  rate 
is  modelled  as  h(t)  =  X(t)  ♦  3  N(t).  We  show  in  Appendix  C  that  the 
characteristic  functional  of  the  hazard  rate  can  be  written  as. 

t 

Ch[ri(-)]  =  exp{  -jX(t’)7l(f)df}  •  (35) 

o 

t  oo  t  +r 

exp{  J2f(t’)dt'li  -  |drb(r)exp[i{^T|(s)ds]  ]  ). 
o  o  t’ 

For  example,  when  the  transit  time  is  exponentially  distributed, 
parameter  p,  the  reliability  can  be  evaluated  as, 

«t)  =  e'xt  expi  -36t/(p*8)  •  3g/(ji*J)2  [l  -  ).  (36) 
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And  when  the  transit  time  is  equal  to  z  with  probability  one,  the 
reliability  is, 

R(t)  =  e"xt  exp{  -tft  -  m  •  (I  -  e'*T)  -  2f(t-r)e"^r},  (37) 

when  t  >  r.  and 

R(t)  =  e'xt  exp{  -2ft  -  2^/3  -  [1  -  e'^) }  .  (38) 

when  t  <  t.  These  formulas  generalize  as  shown  above  for  components 
in  series  and/or  parallel  systems. 

Battle  Damage  Model 

Assume  that  incoming  rounds  hit  a  ship  according  to  a  non-homogeneous 
Poisson  process,  rate  2f(t),  during  a  battle.  The  amount  of  damage  to  the 
ship  per  round  (or,  really,  the  increased  system  stress  resulting  from 
each  hit)  is  given  by  a  random  variable,  3.  with  some  non-negative 
distribution.  Let  the  pdf  of  the  time  to  repair  the  damage  from  each  hit 
be  b(t).  If  the  failure  rate  of  the  total  system  contains  a  term  which 
increases  proportional  to  the  amount  of  un-repaired  damage  (possibly 
because  of  increased  demands  made  on  the  rest  of  the  system),  then  the 
hazard  rate  is. 

n 

h(t)  =  X(t)  ♦  k  I  3j  (0(t  -  tj)  -  0(t  -  tj  -rj)}.  (39) 

j=i 

The  stochastic  portion  of  this  corresponds,  for  example,  to  the  total 
weight  of  customers  in  a  non-homogeneous  M/G/<»  queue,  where  each 
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customer  has  a  random  weight,  3  j.  !n  Appendix  C  we  derive  the 

characteristic  functional  for  this  process.  Applying  that  result  we  have, 
in  general, 

t  t'+r 

R(t)  =  e”At  exp{  -Jdt’2f(r)  [1  -  C^(/i}(s)ds)b(r)dT  (40) 

o  t' 

where  is  the  ordinary  characteristic  function  of  3-  If  ^  is  a  constant 

this  reduces  to  the  previous  model.  Again,  for  m  components  in  series 
replace  X  with  mX  and  evaluate  the  characteristic  functional  for 
Tj (• )  =  m  i.  For  parallel  configurations  Eq.  (C4)  can  be  used  in 
conjunction  with  Eq.  (31). 

One  example  that  can  be  worked  out  in  closed  form  is  when  3  is  an 
exponential  random  variable,  parameter  £,  and  the  system  cannot  be 
repaired  during  the  battle,  i.e.  b(r)  =  5(r  -  <*>).  The  reliability  of  a 
single  component  is  easily  calculated  using  Eq.  (40)  to  be, 


R(t)  =  e"(?*Wt  .  { (e .  t)/e  )£» 
and  for  m  identical  components  in  series  we  obtain. 


(41) 


R(t)  =  .  { (e  ♦  mt)/6  }e*/rn. 


(42) 


Laser-Battle  Model 

We  now  model  the  failures  of  electronic  systems  in  a  laser  battle  as, 

I)  immediate  failures  resulting  from  impulsive  thermal  or  kinetic 
shocks,  or  2)  random  failures  brought  about  by  the  accumulated  thermal 
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stress  from  repeated  hits  by  laser  beams.  The  background  hazard  rate  is 
again  represented  by  \(t).  Assuming  an  exponential  cooling  law,  the 
hazard  rate  resulting  from  laser  hits  at  times  Itj}  is, 

n 

h(t)  =  X(t)  ♦  K<xS(t-t|)  ♦  3  e(t  -  t|)  e"K(t_t j)  },  (43) 

)=> 

where  the  second  term  represents  the  contribution  to  AT  (increase  in 
the  system  temperature)  resulting  from  the  laser  strike  at  time  tj. 

Using  the  same  methods  as  in  the  Appendices  it  can  be  shown  that  the 
characteristic  functional  for  the  stochastic  part  of  this  process  is, 

t  t 

c[t](  )]  =  exp{Jdt'3(t')  I  exp[/ds  t\(s)  r(s  -  t')l  -  U  ).  (44) 

o  V 

where, 

r(s  -  t’)  =  <x  S(s-t’)  ♦  3  e(s  -  t’)  e'K(s'l,).  (45) 

Hence,  the  reliability  is, 

R(t)  =  e"xt  exp{/dt*W)  [  exp[  -<x  -  3/x-(i  -  e"K(t~t>:) )]-!]}.  (46) 

For  m  identical  components  in  series,  seeing  exactly  the  same 
environment,  Eq.  (44)  is  modified  by  multiplying  X,  cx  and  3  by  m.  The 
non-constant  effective  hazard  rate,  h(t),  can  be  read  directly  from 
Eq.  (46),  e.g.  for  m  identical  components  in  series  it  is, 
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h(t)  =  mX(t)  ♦  tf(t)  {  exp[  -m<x  -  m$/x  (1  -  e“K^_t  )  )1  -  1  }.  (46) 

Note  that  %  is  not  multiplied  by  m.  All  components  see  exactly  the  same 
shocks  from  laser  strikes. 

More  generally,  we  could  model  different  types  of  laser  weapons,  ranges, 
atmospheric  propagation,  etc.  by  taking  <x  and  £  as  (correlated)  random 
variables.  The  above  equations  can  also  be  generalized  if  the 
components  are  not  identical  but  have  different  vulnerabilities,  thermal 
conductivities,  etc. 

Vil.  Conclusion 

we  have  presented  a  general  formalism  for  calculating  the  reliability  of 
multi -component  systems  subject  to  a  random  environment.  The 
correlation  in  failures  of  different  components  can  be  accounted  for  by 
calculating  the  system  reliability  in  terms  of  appropriate  combinations 
of  the  characteristic  functional  of  the  random  hazard  rate,  evaluated  for 
■"(•)  =  i*j  {j=l,n;  i  =/(-!)}.  We  have  shown  how  the  method  may  be 
used  in  different  circumstances  by  applying  it  to  a  number  of  different 
models  for  the  random  environment.  Our  results  include  non-constant 
effective  hazard  rates  for  some  of  the  random  environment  models.  The 
intuitively  reasonable  result  that  the  reliability  of  series  (parallel) 
systems  is  greater  (less)  in  a  correlated  environment  than  when  the 
components  see  independent  environments  has  been  demonstrated  for 
some  of  our  models.  Although  we  have  only  considered  series  or  parallel 
systems,  the  application  to  more  complicated  systems  containing 
components  both  in  series  aod  in  parallel  is  straightforward. 
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D.  P.  Gaver.  Jr.  "Random  Hazard  in  Reliability  Problems". 
Technometrics  5.  No.  2,  pp.  21 1-226  (1963).  The  processes 
considered  there  have  independent  increments  (e.g.  the  Gamma 
increments  on  p.  215),  hence  they  yield  a  constant  effective 
hazard  rate  and  exponential  lifetime  distribution. 

2.  A.  Mercer  explores  a  model  of  wear-dependent  failure  rates  that 
leads  to  a  non-constant  hazard  rate  in  "Some  Simple  Wear -dependent 
Renewal  Processes".  J.  Royal  Stat.  Soc.  (B)  23,  pp.  368-376  (1961). 

3.  See  for  example,  R.  E.  Barlow  and  F.  Proschan,  Statistical  Theory  of 
Reliability  and  Life  Testing.  (To  Begin  With  Press.  Silver  Spring.  MD) 
1981. 

4.  M.  S.  Bartlett,  An  Introduction  to  Stochastic  Processes.  2nd  edition, 
(Cambridge  University  Press)  1966. 

5.  D.  R.  Cox  and  V.  Isham,  Point  Processes.  (Cambridge  University  Press) 
1980. 
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6.  E.  Rockower,  "The  Gamma/Exponent ia!  Markov  Process",  NPS  T.R.  1986 

7.  5.  Ross.  Stochastic  Processes.  (Wiley.  New  York)  1983,  p.  2l2ff. 

8.  The  idealized  model  of  a  textile  yarn  discussed  in,  D.  R.  Cox  and  H.  D. 
Miller,  The  Theory  of  Stochastic  Processes.  (Wiley.  New  York)  1965, 
p.  366  ff..  can  be  seen  to  be  equivalent  to  the  M/G/°®  queue. 


Appendix  A.  Derivation  of  the  Characteristic  Functional  for  impulses 
We  now  derive  the  characteristic  functional  for  the  hazard  rate  resulting 
from  random  delta  function  impulses.  Unless  otherwise  noted,  we  define 
the  characteristic  functional  for  the  process  over  the  interval  0  ->  t. 
First,  the  random  part  of  the  hazard  rate,  g(t),  is  defined  more  explicitly 
as, 

n 

g(t')  =  I«j  s(r  -  tj),  (A1) 

where  n  has  a  non-homogeneous  Poisson  distribution,  pn,  with  mean, 
t 

n  =  /  2(t)dt,  (A2) 

o 

the  tj  nave  pdf  given  by  2(0/  n  (as  is  appropriate  for  a 
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non-homogeneous  Poisson  process),  and  o<j  have  an  arbitrary  distribution. 
f(cx j ).  with  ordinary  characteristic  function.  C^.  The  characteristic 
functional  for  g(t')  is, 

t 

Cg[T\(*)l  =  E{  exp(  if'h(t)  g(t)  dt] },  (A3) 

o 

oo  n 

=  I  Pn  ( n  [  JJ at,  «tj)/n  a«|  ««,)  e'  “j  1 ). 
n=0  j=1 

where  the  integration  over  t  j  is  0  ->  t  the  integration  over  <Xj  is  0  ->  <», 

we  have  made  use  of  the  Dirac  delta  function  to  perform  the  integration 
over  r.  The  product  over  j  now  reduces  simply  to  the  expression  in  the 
square  brackets  raised  to  the  nth  power  because  of  the  independence  of 
each  term  in  the  shot  noise- 1  ike  process.  Performing  the  average  over  n 
(yielding  the  standard  result  for  the  generating  function  of  a  Poisson 
distribution)  and  the  integration  over  o<j  (yielding  the  characteristic 

function  for  the  amplitude  distribution)  and  taking  the  limit  T  — >  <», 
results  in. 

t  t 

Ch  th(  )]  =  exp{  iXojT\(t')dt'  ♦  J2f(t’)  (C^  t\(V)  )  -  l]dt’}.  (A4) 

0  0 

where  we  have  included  the  (independent)  characteristic  functional  for 
the  deterministic  portion  of  the  hazard  rate. 
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Appendix  B.  The  Characteristic  Functional  of  the  Gamma  Process6 
Using  the  results  of  Appendix  A,  Eq.  (A4),  the  characteristic  functional 
of  random  delta  function  impulses  (rate  utf(t) )  with  exponential 
amplitudes  (parameter  5).  is, 


00 


Cg[T>( )]  =  exp  {  i  j  dt  otf(t)  ^(t)  /  [S  -  iT|(t)  ] }. 


(Bl) 


-oo 


If  this  shot  noise  is  fed  into  a  first-order  linear  filter  (decay  parameter 
2f(t) ).  the  resulting  process,  x,  satisfies, 

x  ♦  tfx  =  g.  (B2) 


Now,  to  find  the  characteristic  functional  for  the  process,  x(t),  subject 
to  the  initial  condition,  x(0)  =  x0,  use  the  solution  of  the  stochastic 
differential  equation,  (letting  %  be  constant  to  make  it  easier  to  follow 
the  derivation)  in  the  definition, 


Cx[t]OJ  =  E{  exp[  i/r|(t)  x(t)  dtj }.  (B3) 

o 

Or, 

=  exp[ix0Jdt  Ti(t)e‘^1]  •  expii/dt  T^Ue'^/dt'e^  e(t-t’)  g(t’) }, 

where  we  have  used  the  properties  of  the  Heaviside  unit  step  function, 
e(t-t').  so  that  all  the  integrations  are  over  0  ->  «.  thus  making  it 
easier  to  interchange  the  order  of  integration  to  give, 
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Cx[tl(-)1  =  (B4) 

exp[ ix0/ dt  T\(t)e'^t]  •  expfi/dt'  g(t‘)  C/dte"^"*  ^  e(t-t')  T|(t)| }. 

The  second  exponential  term  is  now  of  exactly  the  same  form  as  the 
characteristic  functional  for  g(t),  with  the  expression  in  square  brackets 
replacing  *r\(t’).  Hence,  making  use  of  Eq.  (Bl)  we  have. 

00 

Cx[t](-)1  =  exp[ix0{dt  T|(t)e"^t]  *  (B5) 

o 

00  00  00 

exp  { i / dt '  utf  [{dt  e~^_t  ^  n(t)l  /  [S  -  i  [Jdte"^-t  ^q(t)}  I }, 
o  t'  t* 

where  we  have  used  the  properties  of  the  step  function. 

To  determine  the  marginal  distribution  of  this  process  use 

Tj ( t )  =  Tj o  8(t-t0),  which  recovers  the  ordinary  characteristic  function 

for  x(tg).  A  straightforward  calculation  yields, 

-xft0)‘^0]  =  expliTl0x0e'yto]  •  {  [S  -  i*q0  e"^to]  /  [s  -  i%]  }u.  (B6) 


This  is  the  characteristic  function  of  a  r.v.  that,  with  probability 
c*  -  e"^,  has  the  value  x0e~^  ,  and  with  probability  (l-c*)  is  the  sum  of 
x0e~^  and  a  r.v.  with  exponential  distribution  of  parameter  S. 

Now,  taking  to  -->  ®°  Eq.  (B6)  yields  the  marginal  distribution  of  the 
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steady  state  process, 


Cx(oo)[Tiol  =  (8/  [8  -  i-noW0  .  (B7) 

i.e.  the  characteristic  function  of  a  Gamma  distribution,  as  promised. 

Appendix  C.  The  Characteristic  Functional  of  an  M/G/qq  Queue,  and 
Extensions 

The  M/G/oo  queue  may  be  modeled  as  shot  noise  impulses  of  unit 
magnitude8  (  "customers"  arrive  with  rate  2f(t),  n  is  again  Poisson  ), 

n 

g(t)  =  I  s(t  -  tj).  (CD 

J=' 

which  have  been  filtered  through  a  linear  system  with  a  random  response 
function  given  by. 

r(t  -  tp  =  @(t  -  tj)  -  e(t  -  tj  -  Tj).  (C2) 

i.e.  a  unit  height  pulse  with  random  duration,  tj  (the  random  "service" 
time).  Hence,  for  =  I.  the  number  in  the  system  3t  time  t  is, 


n 

N(t)  =  2  3j  (e(t  -  tj)  -  e(t  -  tj  -tj)}.  (C3) 


we  have  included  the  factor,  so  that  we  can  allow  the  system 


response  to  also  have  a  random  height,  corresponding  to  a  random  weight, 
(or  damage,  system  stress,  pollution,  etc.)  of  each  customer.  Using  the 
definition  of  the  characteristic  functional  and  following  the  same  steps 
as  in  the  previous  Appendices,  we  obtain, 

oo  n  °o 

I  pn  {  n  [  JJJdtj3(tj)/n  d$j  f(0 p  dtjb(rj)exp|i0jJds‘n(s)  ©(s.tj.rp] }. 
n=0  j=!  o 

where  0(s,tj,rj)  =  0(t  -  tp  -  e(t  -  t  j  -tp.  Again  using  the 
independence  of  each  term  in  the  product,  and  the  fact  that  they  are  all 
the  same,  we  can  perform  the  summation  over  pn  (Poisson),  and  the 

average  over  3  (obtaining  the  ordinary  characteristic  function,  ). 
Finally,  the  characteristic  functional  may  be  written, 

t  t'+r 

C[t|(-)I  -  expl -JdtW)  [l  -  C^(jT\(s)ds)b(r)dT  ] }.  (C4) 


o 
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